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The damping of the oscillations of a small permanent magnet (spherical shape, radius 0.1 mm) 
levitating between two parallel YBCO surfaces is measured as a function of oscillation amplitude 
and temperature. The losses in the samples (epitaxial thin films, bulk granular and bulk melt- 
textured) are analyzed in terms of oscillating shielding currents flowing through trapped flux lines 
whose motion gives rise to electric fields. We find dissipation to originate from different mechanisms 
of flux dynamics. At small amplitudes there is a linear regime described by a surface resistance 
varying from 10 -9 fl for the bulk samples down to 10 -13 for the thin films at low temperatures. 
With increasing amplitude various nonlinear regimes are observed, firstly collective pinning with 
diverging energy barriers, secondly in bulk samples above 85 K hysteretic damping, and finally in 
thin films exponentially large losses which can be described by pinning energies vanishing linearly 
at large currents. 

PACS numbers: 74.60.Ge;74.76.Bz;74.80.Bj. 



I. INTRODUCTION 



A permanent magnet levitating above, or suspending 
below, a superconducting surface has an equilibrium po- 
sition which is determined by magnetic forces due to 
screening currents and trapped flux.El Except for a small 
logarithmic long-term relaxation of the magnetization of 
the superconductor this position is stable. The system is 
in a metastable state with stationary flux lines and con- 
sequently without any decay of the dc supercurrents, i.e., 
without any energy dissipation. However, oscillations of 
the levitating magnet about its equilibrium position give 
rise to ac magnetic fields at the surface of the supercon- 
ductor which modulate the shielding currents. Flux lines 
are now moved back and forth periodically as a result 
of the oscillating component of the Lorentz force. This 
leads to energy dissipation which we can determine from 
measurements of the oscillation amplitude of the magnet 
as a function of an external driving force. The measured 
losses can be attributed to various regimes of flux motion 
in the superconductor. 

In our present work we have investigated the energy 
dissipation when the magnet is levitating between two 
horizontal YBa2Cu307_5 (YBCO) surfaces either of bulk 
material (both sintered and melt-textured) and of epitax- 
ial thin films. In all samples and at all temperatures we 
find that at small amplitudes of the ac shielding currents 
(~ 10 A/m) the dissipation can be described by a surface 
resistance which varies from 10 -9 fi for the bulk sam- 
ples down to below 10 -12 fi for the thin films at 4K. At 
larger amplitudes nonlinear behavior is observed which 
can be attributed to thermally activated flux flow with 
pinning energies whose diverging current dependence is 
typical for collective pinning. In addition, for a given 
current, the pinning energies have a linear temperature 
dependence which is an indication of flux creep. Upon 



further increase of the drive and at temperatures above 
85 K we find that the dissipated power in the bulk sam- 
ples varies with the cube of the amplitude. We discuss 
this regime in terms of hysteretic dissipation having a 
coefficient which diverges near T c . In the films, however, 
a quite different behavior is observed instead, namely, an 
exponential increase of the dissipation which ultimately 
exceeds even the bulk values. In this regime we can de- 
scribe the data by a Kim- Anderson type model of flux 
motion with linearly vanishing pinning energies at large 
currents. In summary, by varying the drive of the oscil- 
lating magnet by five orders of magnitude we can study 
different regimes of vortex dynamics in YBCO at a given 
temperature. Some of our preliminary results on granular 
samplesn and on thin filmsa have already been published 
in short form. In this article a more detailed analysis and 
a more complete account of our work is presented. 

This article is organized as follows. In the next Chap- 
ter we describe our experimental method, and in Ch.III 
the model used for data analysis is outlined. In Ch.IV we 
present the results of our experiments for the thin films 
and for the bulk samples. Chapter V concludes our work. 

II. EXPERIMENTAL METHOD 

Our experimental method consists of placing a mag- 
netic microsphere made of SmCos (radius 0.1mm, mass 
to = 5.2 • 10 -8 kg) between two horizontal surfaces (spac- 
ing 1mm, diameter 4 mm) of a parallel-plate capacitor 
made of YBCO. This is an extension o£-|earlier work, 
where the capacitor was made of niobiumQ and the lev- 
itating sphere was used for hydrodynamic studies in su- 
perfluid helium.El In Fig. 1 a schematic of our electronics 
is shown. Before the capacitor is cooled through T c we 
apply a dc voltage of about 800 V to the bottom elec- 
trode. Therefore, the magnet carries an electric charge 
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q of about 2 pC when levitating. The dc voltage is then 
switched off and oscillations of the magnet can be excited 
with an ac voltage U ac ranging typically from 0.1 mV to 
20 V and having a frequency at the resonance of the os- 
cillations (« 300 Hz). These vertical oscillations induce 
a current qv/d in the electrodes, where v is the veloc- 
ity of the magnet, which is measured by an electrometer 
connected to a lock-in amplifier. The capacitively cou- 
pled signal due to the driving voltage U ac which is super- 
imposed upon the signal from the oscillating magnet is 
electronically nulled off resonance by adding the inverted 
pickup of appropriate amplitude. For a given amplitude 
of the driving force F — ql/ ac /d we measure the maxi- 
mum signal when slowly sweeping the frequency towards 
resonance. Because at large oscillation amplitudes the 
restoring forces become nonlinear, the resonance curves 
are usually hysteretic and care has to be taken to achieve 
resonance. We therefore also monitor the phase of the 
signal with respect to that of the drive. In the Appendix 
a short discussion of some of the elastic properties of this 
oscillator is presented. Further data are taken either by 
varying the driving force at constant temperature or by 
sweeping the temperature at a fixed driving force while 
keeping the oscillator at resonance by adjusting the fre- 
quency of the drive. 
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FIG. 1. Schematic of the experimental setup. The elec- 
trically charged magnetic microsphere is levitating inside the 
superconducting YBCO capacitor. 

In order to infer the velocity amplitude from the mea- 
sured current amplitude and the driving force from the 
ac voltage we determine the charge q by recording a reso- 
nance curve at very small drive (which is not hysteretic) 
and by analyzing its width (which is given by the damp- 
ing coefficient 7) and its height which then determines 
the charge. Alternatively, in particular at low damp- 
ing, we measure the time constant r of the freely decay- 
ing oscillations from which we obtain 7 = 2m/ t. Com- 
paring this result with the driven response at the same 
temperature gives the charge. The driving forces range 
from 10~ 13 N to 10 _7 N, the velocity amplitudes from 
0.1 mm/s to 50 mm/s, and the oscillation amplitudes 



a = v/uj from 50 nm to 30 /xm. The dissipated power 
varies from 10 -9 W down to below 10 -17 W. 

The YBCO we used for the superconducting capaci- 
tor was either home-made sintered materials (thickness 
2 mm) or melt-textured sanaplesa of the same geome- 
try, and epitaxial thin filmscl having a thickness of ei- 
ther 450 nm or 190 nm laser deposited on insulating sub- 
strates of SrTiC>3 (lower electrode) and Y203-stabilizcd 
Zr02 (upper electrode) both cut with an inclination of 
2 degrees off the (001) orientation. The inclination im- 
proves both flux pinning and, due to initiating a terrace 
growth, crystallinity of the films. The lower electrode is 
protected by a 20 nm thick PrBa2Cu307 epitaxial layer. 
The YBCO films have an extremely sharp superconduct- 
ing transition at T c — 89.5 K measured by an ac sus- 
ceptibility method, the complete transition width being 
0.1 K. Further details of the thin film production can be 
found in Ref. [?| 



III. DIPOLE MODEL 



The oscillations of the magnet give rise to ac surface 
currents which exert an oscillating Lorentz force on the 
flux lines. This causes energy dissipation which can be 
determined from the damping of the oscillations cither 
by measuring the amplitude as a function of the driving 
force or by recording a free decay. In this section a brief 
discussion of a model used for a quantitative analysis of 
our data is given. 
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FIG. 2. Geometry of the dipole moment above the super- 
conducting surface and its image below. 

In our experiment the superconductors are cooled in 
the field of the magnet and therefore flux is trapped in 
the superconducting electrodes. We estimate the density 
of the trapped flux to be n$o ~ 1 mT corresponding to 
a flux line density n ~ 5 • 10 n m -2 . After cooling we 
find that the magnet always levitates near the middle 
between the superconductors, i.e., at a levitation height 
ho = 0.5mmJ3 In this position the field of the magnet 
at the surface of the superconductors is below B c ± and 
for a quantitative analysis we approximate the surface 
shielding currents induced by the magnet by an image 
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dipole model. In this model the magnetic field on either 
surface is given by the dipolar fields of the magnet and 
the first image dipole (Fig. 2), neglecting images further 
away and also modifications of the field due to trapped 
flux lines. The magnetic moment of the magnet calcu- 
lated from the magnetic remanence of the material and 
the radius is Mq = 4 • 10 -6 Am 2 . It is simple to calcu- 
late the dc sheet current on the surface for a particular 
orientation of the dipole with respect to the surface.EJ 

The dc current distribution J = {J x {x, y, h), J y (x, y, h)) 
(x and y are coordinates on the surface) is given by 
J x = —2H V and J v — 2H X with 



H x = - — - [xy cos a — hy sin a) , 
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where r 2 = x 2 + y 2 + h 2 . The parallel orientation of 
the dipole (a — 0) has the lowest energy but because 
of trapped flux we cannot rule out that some other ori- 
entation applies. Although the variation of the surface 
current is quite complex and changes drastically with the 
orientation of the dipoletL9, the maximum value J max is 
fairly insensitive to the orientation: J max = 4400 A/m 
when the dipole is perpendicular and J max = 5100 A/m 
for the parallel orientation at h = ho. In either case 
the magnetic field is below H c i, but some flux remains 
trapped because of the field-cooled situation. 

The ac sheet current A J due to the oscillations of the 
magnet about its position at ho is 
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acos(uit), (3) 



where a is the oscillation amplitude. Again the cur- 
rent distribution (shown in Fig. 3) is complicated but 
the maximum amplitude of AJ is insensitive to the ori- 
entation of the dipole. For an oscillation amplitude 
a = 0.5 fim corresponding to a velocity amplitude of 
v rs 1 mm/s (at 300 Hz oscillation frequency) a maximum 
current amplitude of 14 A/m is found for the perpendic- 
ular orientation and 15 A/m for the parallel one. The 
current is localized in a small area of the superconduct- 
ing surfaces below and above the sphere, i.e., dissipation 
occurs only here. 

This ac sheet current is actually an oscillating current 
density j in the samples. Flux lines will be set into mo- 
tion thereby creating an electric field E = Vf-n^o, where 
vj is the flux line velocity. Because j and E are par- 
allel the dissipation per unit volume is given by j ■ E. 
In bulk samples the shielding current is flowing within 
the London penetration depth where it decays exponen- 
tially, whereas in thin films the penetration depth can 
exceed the sample thickness. (Note that the Campbell 
penetration depth is an order of magnitude smaller be- 
cause of the low trapped flux density and therefore can 
be neglected.) For simplicity, we use an average current 



density (j) = AJ/5 inside the samples, where 6 is either 
the penetration depth in case of bulk samples or the film 
thickness. We may then describe the dissipation per unit 
area by AJ • E, where E(j) depends on (j). 
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FIG. 3. Vector plots of the ac surface currents calculated 
from Eq. ^[ a) for parallel and b) for perpendicular orienta- 
tion of the dipole. In our case h — 0.5 mm and therefore the 
displayed surface area is 0.5 x 0.5 mm 2 each. 



We calculate the stationary velocity amplitude of the 
magnet v(F) = wa(F) by employing the energy balance 
between the gain from the driving force (which at reso- 
nance is in phase with the velocity) and the loss per cycle 
(period r) on both surfaces: 
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Because A J oc a it is evident from Eq. [| that a linear de- 
pendence of E(j) corresponds to a linear response v(F) 
while a nonlinear response implies a nonlinear current- 
voltage relation. Secondly, it is obvious that the driving 
force is a measure of the electric field, i.e., of the vortex 
velocity. In the following chapter we present our experi- 
mental results and analyze them on the basis of Eq. 0. 



IV. RESULTS AND DISCUSSION 

We first describe and discuss the data obtained with 
the thin film samples which are very similar for both the 
190 nm and the 450 nm thick films. Then, the measure- 
ments on the bulk samples are presented which are very 
similar for both granular and melt-textured YBCO. 



A. Thin films 



A summary of our data obtained with the 450 nm films 
above 77 K is shown in Fig. 4. At small driving forces we 
find a steep initial increase of the velocity amplitude of 
the magnet, whereas at driving forces above ca. 10 nN 
there is only a logarithmic dependence. This implies an 
exponential increase of the dissipation. At lower temper- 
atures the velocity amplitudes increase very rapidly with 
drive as shown in Fig. 5 with a largely expanded abscissa. 
Further increase of the driving force amplitude gave un- 
stable oscillations. Only at very small amplitudes we 
do find a linear regime between the driving force ampli- 
tude F and the velocity amplitude v, i.e., F = 71; with 
a strongly temperature dependent coefficient j(T), see 
Fig. 6. While 7 decreases considerably below the critical 
temperature T c , it obviously levels off at 4K. Here the 
Q-factor 77107/7 which-we attribute exclusively to losses 
in the superconductom exceeds 10 6 . The steep drop be- 
low T c follows approximately a (T c — T)~ 2 law as can be 
seen in the insert of Fig. 6. 

In the linear regime F = 77; the induced electric field 
in the superconductors is proportional to the surface cur- 
rent amplitude, E = R S AJ, and the loss (r.h.s. of Eq. |^) 
can be described by a surface resistance R s oc 7, namely 



R, = 



where 



5 = 2 




dxdy. 



(5) 



(6) 



Here we assume that R s is a constant on the surface of 
both superconductors. The surface integral S is evalu- 
ated numerically. In Fig. 6 the R s values are plotted on. 
the right scale for the parallel orientation of the dipole.EJ 
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FIG. 4. Resonant velocity amplitude of the oscillating 
magnet as a function of the driving force at various temper- 
atures. Note the steep increase at small driving forces and a 
logarithmic dependence at large drives. The lines are fits of 
our model to the data, see text. An almost linear dependence 
observed with melt-textured bulk samples (from Fig. 11) is 
shown for comparison. 
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FIG. 5. Velocity amplitude at very small driving forces in- 
dicating the linear initial increase at various temperatures. 
The lines are fits of our model to the data, see text. 

From the surface resistance we can infer the imaginary 
part of the penetration depth of the ac field being given 
by A" = R s /hogj, In our case (uj/2ir w 300 Hz) we find 
A" to drop from 84 nm at 85 K down to only 0.2 nm 
at 4 K. A steep decrease of the surface resistance with 
temperatuKe,has been observed at radio and microwave 
frequencies^ 2 ] and is discussed in terms of various mod- 
els of flux dynamics. Ln_the limit of very low frequencies 
it has been calculatedEj that the complex resistivity p v 
due to vortex motion is greatly altered when flux creep 
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is included. Instead of a purely imaginary resistivity the 
result has a real part given by Re(p„) = pt ■ e where 
Pf = B^q/t] is the usual resistivity for free vortex flow 
(rj is the vortex viscosity) and e{v) = Iq 2 (v) <C 1, where 
Io is a modified Bessel function and v = U/2kT depends, 
on the pinning energy U (T) . There is ample evidenceo 
that below 70 K v is a temperature independent quantity 
ranging from 3 to 9. In fact, for v = 9.5 our results would 
agree very well with R s — Re(p v )/t (t is the film thick- 
ness) if a trapped flux density of 1 mT and typical values 
for r)(T) are assumed, see Ref. [l4] for a compilation of 
viscosity data. This value for v appears quite reasonable 
because, as will be shown below, in the nonlinear regime, 
where v begins to drop with increasing current, we find 
v = 8 at a current amplitude of 50A/m. It should be 
mentioned that recent work oiuthe surface impedance of 
conventional superconductors^ has shown that a more 
complicated electrodynamical description is required. If 
this theory were directly applied to our experiment (thin 
films, low frequency, low vortex density) one would ob- 
tain the simple relation R s oc (T) which would agree 
both in temperature dependence (including the region 
near T c ) and in order of magnitude with our data. How- 
ever, this remains to be justified and therefore at present 
is only a speculation. 
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FIG. 6. Linear coefficient 7 = F/v from Fig. 5 and surface 
resistance R s from Eq. [|. Open symbols are data obtained 
with the thinner film 5 = 190nm. The solid line is the same 
as in the insert where the data are plotted vs T c — T. 

In Fig. 7 the data of Fig. 4 are plotted on a logarith- 
mic scale for the driving force F. Above 10 nN we find 
a logarithmic dependence v oc \ii(F/Fq) which implies 
an exponential increase of the losses. For driving forces 
between 0.1 nN and 10 nN the amplitudes appear to be 
reduced indicating an additional loss mechanism in this 
regime. Furthermore, it is noteworthy that obviously the 
driving force and hence the electric field, or vortex veloc- 
ity, determines which mechanism dominates and not the 
oscillation amplitude, or shielding current. 

In the nonlinear regime of v(F) the electric field in 



Eq. ^ must have a nonlinear current dependence. This 
is usually described by thermally activated flux motion 
E(j) oc exp(— U(j)/kT) with current dependent pinning 
energies U(j), for a review see Ref. [r| Essentially two 
different dependences are being considered, firstly barri- 
ers diverging for j — ► 0: 

U(j) = U (jY, (7) 

where < p < 1. Vortex glass and collective creep mod- 
els yield this behavior which implies zero dissipation in 
the limit j — + 0. Secondly, for large currents, one con- 
siders barriers vanishing at a critical current density j c 
as 

U(j) = U (l-j-), (8) 

which leads to an exponential dissipation for large cur- 
rents (Kim-Anderson model), E oc sinh(j/j'i) with ji = 
JckT/Uo. 
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FIG. 7. Data of Fig. 4 plotted on a logarithmic scale for 
the driving force. The solid lines are fits of Eq. ^ to the data, 
see text. For the 81K curve the contributions of three differ- 
ent mechanisms of dissipation are indicated separately: 1 is 
the linear one (which is an exponential curve because of the 
logarithmic abscissa), 2 is the contribution from Eq. |^, and 3 
is from Eq. ||. 

We fit the complete v(F) curves by including both non- 
linear mechanisms (Eqs. |7] and ^|) and a linear term ad- 
ditively in the energy balance, see Eq. [|. An example 
is shown in Fig. 7, where the separate contributions and 
their sum are compared with an experimental curve. The 
linear term (surface resistance) contributes only at very 
low amplitudes. In the intermediate range the glass (or 
creep) term (Eq. |t] with // = 0.2) is dominant. Only at 
the largest drives the Kim-Anderson term (Eq. ^) be- 
comes relevant. It is this mechanism which leads to the 
observed slow logarithmic increase of the oscillation am- 
plitudes. From these fits we can determine the barriers 
U(j) at various temperatures, see Fig. 8. It is a peculiar- 
ity of Eq. ||that only the constant slope — dU /dj = Uo/j c 
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can be obtained reliably from a fit to the data because 
the quantity Uq determines only the prefactor which we 
do not evaluate for it contains other unknown quantities, 
e.g. the density of the trapped flux n or the prefactor 
of the vortex velocity v / . From Fig. 8 it is evident that 
at large currents and high temperatures the barriers of 
Eq. H ultimately drop below those of Eq. [7| For the ex- 
ponent [i determining the barriers we find fi = 0.2 which 
is close to the theoretical value 1/7 for single vortex pin- 
ning which is reasonable for the small vortex densities in 
our experiment .E3 We should mention that the exponen- 
tial fitting parameters are not affected by the particular 
current distribution used (parallel or perpendicular orien- 
tation of the dipole), only prefactors are different because 
the parallel dipole leads to a lower total dissipation than 
the perpendicular one. 
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FIG. 8. Dependence of the pinning potential U on the av- 
erage current density in the film obtained from the fits of 
our model to the data: curved lines are Eq. [?], straight lines 
are Eq. |^ at the following temperatures (from top to bottom): 
78 K, 81 K, 83 K, 84.5 K, 85.5 K. Note that from Eq. | we can 
determine only the slopes Uo/jc', in the figure we arbitrarily 
set U = 4 • 10 3 K 

It is interesting to note that all three mechanisms of 
dissipation are found to coexist at the same temperature 
but dominate in different regimes of the driving force. 
Furthermore, we find that the current j\ which sets the 
scale for the exponential dissipation extrapolates to zero 
at 86.5 K (well below T c = 89.5 K) which might indi- 
cate an onset of free flux flow where the magnet loses its 
lateral stability. 

Towards lower temperatures only the linear term and 
the glass term contribute as no data could be obtained 
at larger drives because the oscillations became unstable. 
For temperatures below 60 K we find that the evaluated 
barriers decrease nearly linearly with temperature, see 
Fig. 9. We find U = 16T and U = 10T (i.e., v = 8 and 
5) for current amplitudes of 50A/m and 500 A/m, re- 
spectively. The ratio 1.6 is given by the factor 10 ' 2 , see 
Eq. |7|. The linearity of U (T) is again a signature of flux 
creep, but now with a current dependent creep factor. 



The values for v are determined by the logarithm of var- 
ious properties of the pinning centers and of the inverse 
frequency. For this reason our values obtained at 300 Ha 
are larger than those obtained at microwave frequenciesEj 
but lower than those obtained from magnetic relaxation 
measurements. Approaching T c the effective barrier fi- 
nally decreases (see Fig. 8) because the quantity jo p. 
Eq. [?] which is related to the critical current density££l 
goes to zero. 
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FIG. 9. Temperature dependence of the effective barriers 
(Eq. for two different current amplitudes. Below 60 K the 
dependence is linear for both currents with slopes 10 and 16, 
respectively. The solid lines are merely a guide to the eye. 



B. Bulk samples 

We will now give an overview of the measurements with 
YBCO bulk samples and compare these results with the 
thin films. We have found that the dissipative behav- 
ior of melt-textured and sintered bulk samples is quite 
similar. Measurements on the sintered samples have al- 
ready been publishedtl and therefore we will focus on our 
melt-textured results. We emphasize, however, that our 
present interpretation of the data is more complete and 
to some extent makes our earlier analysis obsolete. 

Figures 10 and 11 show the velocity amplitude of the 
magnet as a function of driving force at different tem- 
peratures for melt-textured and sintered samples. At 
4 K we find a nonlinear increase of the amplitude even 
at small amplitudes followed by a nearly linear depen- 
dence for large drives. This initial nonlinear increase is 
less pronounced at higher temperatures where we find 
approximately linear curves up to 85 K. For higher tem- 
peratures additional dissipation arises and above 88 K 
the dependence is simply v oc V~F indicating a quadratic 
damping, i.e., losses proportional to the cube of the am- 
plitude. Because of the various mechanisms of dissipation 
(to be discussed below) dominating at diffent tempera- 
tures and amplitudes, there exists a temperature range 
below T c where the oscillation amplitudes surprisingly 
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increase with temperature at constant drive, correspond- 
ing to a minimum of dissipation. This can be seen more 
clearly in Fig. 12 where the velocity amplitudes are shown 
as a function of temperature for various constant drives. 
We note that for sintered YBCO the maximum is much 
more pronounced, see Ref. g. 
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FIG. 10. Velocity amplitudes as a function of the driving 
force for sintered and melt-textured samples at 4K and 77 K. 
The lines are fits of our model to the data. 
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FIG. 11. Velocity amplitudes versus driving force at tem- 
peratures above 85 K where damping becomes quadratic lead- 
ing to a square-root dependence. The lines are fits of our 
model to the data. 

Considering the initial nonlinear increase of the am- 
plitude at 4K followed by a nearly linear rise for larger 
drives we find that this behavior can be well fitted by 
adding a linear dissipation and the collective flux creep 
(Eq. 0), using an exponent fj, = 0.14. Because Eq. 
describes the data over almost the entire range the ex- 
ponent fi can be determined with much better accuracy 
than with the thin-film samples: increasing \x to 0.2, e.g., 
leads to a clearly poorer fit. From the linear term (domi- 
nating only at the smallest amplitudes) we infer a surface 
resistance which is of the order of R s ~ 10~ 10 SI at 4K 



rising to 10 -8 ft at 85 K. Comparing these results with 
our data on YBCO thin films we note that at 78 K R s 
has increased by three orders of magnitude. 

At 4K the barrier in the nonlinear regime is much 
smaller than at high temperatures (there are no experi- 
mental data between 4K and 77 K). Above 77 K the bar- 
riers decrease with increasing temperature. These results 
are similar to those obtained with the films except for the 
fact that the barriers in the bulk are a factor of 2 smaller 
than in the films. 
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FIG. 12. Temperature dependence of the velocity ampli- 
tude at various constant driving forces. Note the maximum 
at about 88 K (i.e., minimum of dissipation). 

For temperatures above 85 K quadratic damping 
rapidly increases and eventually becomes to be the only 
relevant dissipation above 88 K. We attribute this to a 
mechanism proposed by MelvilleHl for conventional su- 
perconductors in which single vortices are considered be- 
ing trapped perpendicular to the surface of the super- 
conductor. When an external field smaller than H c i is 
applied shielding currents are induced at the surface pro- 
ducing a Lorentz force on the trapped vortices. If this 
force is strong enough to depin the vortices they will 
bend in the direction of the applied force, while they 
stay pinned in the interior of the bulk where no current 
is flowing. The amount of bending is given by the equi- 
librium between the Lorentz force and the elastic force 
of the vortex. When the applied field is reduced and fi- 
nally reversed the vortex moves back in an irreversible 
way, so that the flux component parallel to the surface 
passes through a hysteresis. This leads to dissipation 
where the loss per cycle and surface area is given by 
B H 3 1 H c i j c , where H is the amplitude of the external 
field and j c the depinning current density. In our exper- 
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iment, H oc a and therefore, with Eq.|f|, this mechanism 
leads to quadratic damping F oc v 2 having a coefficient 
which diverges at T c as was found experimentally already 
for sintered samples a 

It should be noted that the well known Bean critical 
state also leads to quadratic damping and is used suc- 
cessfully in the interpretation of most other levitation 
experiments involving larger magnets. However, in our 
experiment the magnetic field is low and a critical state 
is not reached. Also, in the model of Melville the dis- 
sipation is proportional to the density of trapped flux 
B which makes it plausible that the losses in the melt- 
textured samples are larger than those in the sintered 
samples. 

Besides the much larger damping found in bulk sam- 
ples, there are two main features distinguishing the bulk 
results from the thin film ones: Firstly, for bulk samples 
the damping becomes quadratic approaching T c so that 
the amplitudes exceed those for thin films (Fig. 4) where 
the dissipation grows exponentially with amplitude. We 
believe that the reason for this is the sample thickness. 
In bulk samples the vortices stay pinned in the interior 
and move just within the penetration depth. In thin films 
the vortices are exposed to the driving current over their 
entire length. This leads to high flow rates of the vortices 
in the films whereas in bulk samples the flow rates are 
limited due to pinning deeper inside the samples. 

Secondly, for bulk samples there is a minimum in dis- 
sipation below T c while for thin films the damping grows 
monotonically with temperature. This minimum might 
be due to a softening of the flux line lattice leading to 
improved pinning because the vortices can adjust more 
effectively to a state of lower total energy. At still higher 
temperatures pinning becomes less effective, hysteretic 
damping occurs and dissipation diverges approaching T c . 
In thin films the exponential growth of dissipation dom- 
inates and no minimum in damping occurs. 



V. CONCLUSION 

In summary, the damping of the oscillations of the lev- 
itating magnet is due to energy losses caused by flux mo- 
tion of various dynamic regimes. At sufficiently small am- 
plitudes linear friction is observed which can be described 
by a surface impedance having particularly low values for 
thin films at low temperatures. The temperature depen- 
dence of the surface impedance is discussed in terms of 
a vortex resistivity caused by flux creep. Towards larger 
amplitudes we find losses to originate from the motion of 
collectively pinned vortices. From our analysis pinning 
barriers are obtained having the typical inverse current 
dependence. The effective barriers increase linearly with 
temperature up to 60 K as expected from creep effects. 
Finally, at even larger amplitudes, the losses in the films 
grow exponentially and eventually exceed those in the 
bulk samples. This seems plausible because in the films 



the flux lines are exposed to large current densities along 
their entire length, whereas in the bulk the Lorentz forces 
decay with increasing distance from the surface. 

When comparing our experimental method (which was 
motivated by an interest in the physics of superconduct- 
ing levitation) with more standard techniques for inves- 
tigating vortex dynamics in superconductors, like mea- 
surements of the ac susceptibility, of the current volt- 
age characteristics ,of the magnetic relaxation, or me- 
chanical oscillators containing superconducting samples, 
we note several obvious differences some of which are 
disadvantageous while others are of advantage. Firstly, 
our method is not affected by demagnetization factors or 
edge effects, the geometry is just a magnetic dipole in 
the middle between essentially infinite and flat supercon- 
ducting surfaces. Secondly, it is a resonance method with 
high sensitivity and wide dynamic range: The dissipated 
power varies from 1CP 9 W down to 10 -17 W and the elec- 
tric fields caused by vortex motion range from I0~ 4 V/m 
down to 10 -11 V/m. Because of the resonance method, 
though, the frequency cannot be varied. The magnetic 
field is inhomogeneous and rather small, and the distri- 
bution of trapped flux lines is unknown. Nevertheless, 
the damping of the oscillating magnet can be described 
quantitatively by various mechanisms of flux dynamics 
with flux creep being an important issue. 
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APPENDIX: ELASTIC FORCES ON THE 
OSCILLATING MICROSPHERE 

The elastic forces on the levitating microsphere were 
investigated by measuring the amplitude dependence of 
the resonance frequency co(a) and also of the second and 
third harmonic amplitude. In addition, the effect of a 
large dc voltage (up to ±1000 V) across the capacitor 
on the above quantities was determined. By ancULarge, 
the oscillator displays typical nonlinear properties S3 The 
"backbone curve" u)(a), e.g., is the usual parabola for the 
thin-film samples below 81 K. However, at higher tem- 
peratures and for all bulk samples in particular, stronger 
dependences are found which indicate that a third order 
Taylor expansion of the return forces at the equilibrium 
position is no longer adequate. On the other hand, the 
coefficients of the second and third harmonic arc found 
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to be independent of temperature. Application of the dc 
voltage changes the equilibrium position by few tens of 
micrometers as can be estimated from the stiffness. Fur- 
thermore, the resonance frequency changes because of 
the changing tension, in analogy with a vibrating string: 
u) 2 is proportional to the static force as observed earlier 
with a niobium capacitor .□ The coefficient of the second 
harmonic is strongly affected and can be reduced to zero 
at a particular voltage, which implies a symmetric return 
force at the new equilibrium position. 

These elastic properties may be discussed in terms of 
the image-dipole model presented in Ch. III. Neglecting 
higher order images the resonance frequency is easily cal- 
culated to be bj/2n w 200 Hz if the equilibrium position 
is in the middle between the two surfaces. Gravitation 
and especially trapped flux lead to unsymmetric forces 
and therefore to an unsymmetric position and hence to a 
second harmonic amplitude. The force between a single 
vortex and a magnetic dipole perpendicular above a su- 
perconducting thin film has been considered in Ref. |2(]. 
Because we have no quantitative information on the dis- 
tribution and the orientation of the trapped vortex lines 
in the samples we cannot calculate their contribution to 
the elastic forces on the levitating magnet. Therefore, 
the image-dipole model remains a rather crude descrip- 
tion of the elastic forces. For lateral stability of the mag- 
net, which is a prerequisite of our method, trapped flux 
is clearly necessary. 
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